Energetics and dynamics of H 2 adsorbed in a nanoporous material at low temperature 
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Molecular hydrogen adsorption in a nanoporous metal organic framework structure (MOF-74) 
was studied via van der Waals density-functional calculations. The primary and secondary binding 
sites for H2 were confirmed. The low-lying rotational and translational energy levels were calculated, 
based on the orientation and position dependent potential energy surface at the two binding sites. A 
consistent picture is obtained between the calculated rotational-translational transitions for different 
H2 loadings and those measured by inelastic neutron scattering exciting the singlet to triplet (para 
to ortho) transition in H2. The H2 binding energy after zero point energy correction due to the 
rotational and translational motions is predicted to be ~100 meV in good agreement with the 
experimental value of ~90 meV. 
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The adsorption of molecules within the nanopores of 
a sparse material and their low-lying excitations provide 
rich phenomena of fundamental interest that are seldom 
explored by first principles methods, even for molecules 
as simple as H2. For example, one can ask how the well- 
known para— ortho transition of H2 survives in a recogniz- 
able form in the presence of the rotational barriers and 
hindrances provided by the adsorbing material, and how 
it changes for different concentrations of H2. The neces- 
sity of proper treatment of van der Waals interactions has 
foiled traditional density functional treatments. Here, by 
using the van der Waals density functional (vdW-DF) of 
Dion et al. [![, we are able to predict a picture of the 
low-lying excitations which provides a credible match for 
the results of inelastic neutron diffraction (INS) for 
different hydrogen loadings in a nanoporous material of 
a type thought to provide a possible path toward future 
hydrogen storage technology. 

The material studied here is a metal-organic frame- 
work (MOF) compound, namely the material that has 
been called MOF-74 0, chosen because of recent in- 
elastic neutron scattering measurements Q with varying 
amounts of H2 adsorption. More generally MOFs are 
a relatively new class of materials which are composed 
of metal clusters connected by organic ligands |4] . The 
search for a MOF structure with high binding strength to 
H2 has become very active recently [5|, |6| . The dynamical 
properties of the adsorbed dihydrogen such as vibrational 
and hindered rotational motion play a significant role in 
determining the H2 binding energy, especially for caged 
structures with nanopores. This area is even less ex- 
plored than the modeling of hydrogen uptake based only 
on the depth of the potential well. An attempt along this 
line studied the rotational transition of H2 adsorbed over 
benzene molecule which is often a fraction of the organic 
linker in MOF materials 0- However, the interaction be- 
tween full MOF and H2 usually differs significantly from 
that between a fraction of MOF and H2 |8|. A very re- 
cent paper studied the rotational transition in another 
MOF (HKUST-1) with generalized gradient approxima- 



tion (GGA) calculations [9(. However, GGA is known to 
fail for cases where the binding is dominated by van der 
Waals interactions, and only a preliminary theoretical 
picture was obtained. 

Here using vdW-DF fij, implemented via the efficient 
algorithm of two of us [IJJ] , we calculated the dynamical 
properties as well as the binding energy of H2 adsorbed 
in MOF-74 structure, which has been demonstrated to 
increase both the H2 affinity [|[ and hydrogen density Q 
due to the unsaturated metal sites 0] . The primary bind- 
ing sites associated with metal atoms are confirmed by 
our calculations with a binding energy of 130 meV, while 
the LDA and GGA results are 230 and 46 meV respec- 
tively Q. The center-of-mass translation and molecular 
rotation of H2 arc found to attribute zero point energy 
(ZPE) of ^20 meV and ~10 meV, respectively, to the 
total energy of the MOF-H2 system. As a result the ef- 
fective binding is about 100 meV, which is close to the 
experimentally measured value of ~90 meV [1, H|. 

Furthermore, we find that the adsorbed H2 is a three- 
dimensional quantum rotor at both the primary and sec- 
ondary binding sites. The original J = 1 triplet state 
splits into three nondegenerate levels, which imply three 
distinct para to ortho excitation energies. A good match 
between the calculated rotational transitions and the in- 
elastic neutron scattering spectra [1] was obtained. As 
suggested in Ref. [TJ, MOF materials may be used to 
create artificial hydrogen nanocages. In addition to the 
confinement and short intermolecular distance, the tem- 
perature dependent preferential adsorbance [12| of ortho- 
H2 due to the splitting of the triplet could give rise to 
unexpected properties as well. 

To locate the binding site near the Zn atom, we used 
the experimental MOF structure. With one H2 per prim- 
itive cell (1/6 H 2 per Zn), we calculated the MOF-H 2 
interaction energy for various H2 positions near the Zn 
atom. The contour map of the interaction energy is 
shown in Fig. [TJ The minimum position has a binding 
energy of 130 meV and its distance to the Zn atom is 
2.9 A while the experimental value is 2.6 A The 
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FIG. 1: Potential energy surface at the Zn site (meV). Top 
panel: cell with energy map plus the location of primary and 
secondary sites at an equivalent place in the cell; the two sites 
have different depths. Bottom panel: blowup of energy map 
in a different orientation. 



over-estimation of the bond length is typical in vdW-DF 
calculations [HI]. 

Along with the strongest binding site near the metal 
atom (hereinafter referred to as the Zn site), a secondary 
site above the oxygen triangle was observed by neutron 
scattering measurements Q ■ We again took the structure 
of Ref . d with a loading of 12 D2/ cell. Instead of sampling 
the potential energy surface for the two H2 manually, we 
performed a relaxation of the adsorbed dihydrogens with 
a fixed MOF host with the force converged to 12 me v/A. 
We find that a H2 indeed positions itself above the oxygen 
triangle when its neighboring metal sites are occupied. In 
other words, a potential minimum is created at the O- 
triangle site due to the adsorbance of H 2 at the metal 
site. The closest distance between sorbate and sorbent 
is 3.4 A, while the distance between the O and Zn sites 
is 3.2 A. These are each 0.3 A larger than the respective 
experimental distances Q. 

With the two binding sites located, we now turn to the 
dynamical properties of the adsorbed H2 at these sites, 
which are important in determining the zero-point en- 



ergy correction to the total binding strength. In general, 
H 2 has six degrees of freedom. We here focus on the 
translational and rotational motions where the H2 bond 
length is fixed at its equilibrium value. Even with fixed 
bond length, the full quantum mechanical treatment of 
the remaining five degrees of freedom is still beyond the 
scope of current letter. We make the approximation that 
rotation-translation coupling may be neglected to lowest 
order, so that the two motions can be treated separately. 
In the absence of intersite dynamical coupling, this would 
be exact to the extent that the H2 motions were small 
enough to be described by the harmonic approximation. 
Neither of these is exact, so that small deviations from 
the rotational spectra we calculate are to be expected, 
even though we use methods that do not explicitly as- 
sume the harmonic approximation's validity. 

Fig.[5]shows the orientational dependence of the inter- 
action energy at Zn site. The maximum point in this 
figure corresponds to H2 bond aligned toward the Zn 
atom, which is the weakest binding orientation, while the 
preferential bond direction of the adsorbed dihydrogen is 
in the plane perpendicular to the Zn-H axis. The O 
site, however, is quite different from the Zn site. The H2 
tends to line up toward the oxygen triangle and the pref- 
erential orientation is inside a conical region (figure not 
shown). The maximum variation of the rotational poten- 
tial at each of the two sites is ^27 meV and ~18 meV, 
respectively. Due to this variation, the adsorbed hydro- 
gen molecule exhibits hindered rotation behavior at these 
two sites. As a result, the energy levels of the free H2 ro- 
tor (E = BJ(J + 1),B = 7.35 meV) will be lifted and, 
in particular, three nondegenerate states emerge from 
what was the J = 1 level with m = 0, ±1, respectively. 
The transition between the singlet and triplet states is 
the well-known para-ortho transition. A spin flip is re- 
quired to satisfy the antisymmetry requirement of the 
total wave function. To compare with INS data 0, we 
calculated the rotational spectrum by solving the rigid 
rotor Schrodinger equation for H2 . We first fitted the ori- 
entational potential E{0, cj)) to a summation of spherical 
harmonics and then diagonalized the Schrodinger equa- 
tion with a spherical harmonics basis. The results are 
shown in Table U where three spectra with different load- 
ings are calculated for Zn site and one set for the O site, 
which is mostly unoccupied at low loading. 

In Table [J we first observe that the ground state ener- 
gies (zero point energies) are about 10 meV at both sites, 
and it increases by about 3 meV at Zn site as the loading 
is increased. The excited state of J — 1 splits into three 
nondegenerate states for all four cases considered theo- 
retically, showing the hindered three-dimensional quan- 
tum rotor character of the adsorbed dihydrogen at the 
two sites. For the Zn site with 1 H^/cell loading, the 
three para-ortho transitions are 10.6, 14 and 21.6 meV, 
which we assign to the main INS peaks at 8.3, 11.1, and 
20.9 meV in Fig. S2 of Ref. 0, whose points we replotted 
in Fig. [3] for easy comparison. In addition to the three 
major peaks, two small peaks at 14.2 and 18.6 meV are 
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FIG. 2: Orientational dependence of the potential at the Zn 
site. The spherical coordinates are defined with respect to the 
axes of Fig. [1] The hydrogen loading is 1 bb/cell. 

also observed experimentally in the INS spectrum. Their 
strength increases substantially as the loading increases 
so that the second site begins to be occupied, and they 
move to slightly higher energies (14.5 and 19.2 meV, re- 
spectively). This result is in surprisingly good agreement 
with our calculations which show that the O site gives 
two transitions at 15.6 and 18.5 meV. Moreover, the ex- 
perimental peak at 11.1 meV significantly increases in 
strength with a shoulder appearing under high H2 load- 
ing; this feature probably consists of two peaks with one 
of them due to the second site and the other due to the 
Zn site. This result is again consistent with our find- 
ings that H2 at the O site has a para-ortho transition at 
~11 meV. From Table Q] we also see that the predicted 
transition energies at Zn site shift as loading increases, 
especially for the Zn-3 transition, which changes from 
^22 to ^26 meV while the other two shift downward 
only by ~1 meV. Experimentally, there is a smaller cor- 
responding shift to about 23 meV for higher loading (9.6 
H2/cell). For the loadings of 7 and 9.6 Eb/cell, we assign 
Zn-3' to the contribution from the sites where the neigh- 
boring O site is empty, while assigning Zn-3 to the case 
where it is occupied. Theoretically there are no empty O 
sites at the 12 H^/cell loading, while presumably at the 
7 H2/cell loading, none of the 5 Zn-3 sites with empty O 
sites would show the large shift that would be caused by 
a filled adjacent O site. 



In Table HT1 we show, for different sites and loadings, the 
translational energies along normal coordinate directions, 
obtained from the forces induced by small displacements. 
The energies were obtained by using these harmonic nor- 
mal coordinates to extend the results approximately into 
the anharmonic quantum-mechanical regime, by first dis- 
placing the H2 along these normal directions and then 
substituting the resulting potential into the Schrodinger 
equation to get the energies Q, and hence the transition 
frequencies. The frequencies determined by this method 
deviate as much as ~10% from the harmonic frequen- 
cies, although for most cases the deviation is less than 



TABLE I: Theory vs. experiment (Fig. [3} for the positions 
(meV) of rotational energy-loss peaks along with the calcu- 
lated rotational zero point energy. (RZPE). 



# Ha/cell 


Method 


RZPE 


Zn-1 


Zn-2 


Zn-3 


1 


theory 


10.1 


10.6 


14.0 


21.6 


4.8 


experiment 




8.3 


11.1 


20.9 


7 


theory 


13.6 


9.2 


13.5 


26.2 


9.6 


experiment 




8.6 


10.6 


23.0 


12 


theory 


12.9 


9.8 


12.8 


26.4 








0-1 


0-2 


0-3 


9.6 


experiment 




11.2 


14.5 


19.2 


12 


theory 


9.6 


11.0 


15.6 


18.5 



5%. Table HIl shows that the translational frequencies are 
remarkably close to the rotational transitions shown in 
Table U and could lead to a complex INS spectra. How- 
ever, the scattering probability for para-ortho transition 
is proportional to the incoherent cross section (<7j) of 
atomic hydrogen while the para-para scatte ring is pro- 
portional to the coherent cross section (cr c ) [lj|. These 
cross sections have been tabulated [l5[ and one finds 
that Oi > 40cr c . Therefore the INS spectra would be 
dominated by the para-ortho transitions with little di- 
rect contribution from translational motions. Simultane- 
ous excitation of rotational and translational excitations 
may appear as shoulders or sidebands of the major para- 
ortho transition peaks, which are not necessarily weak. 
Although two possibilities are labeled (T-l and T-2) in 
Fig. [3l an unambiguous assignment is difficult. The ma- 
jority of such double excitations would be predicted to 
occur at frequencies higher than those measured, or else 
hidden by the strong Zn-3 structures assigned to rota- 
tional excitations. Spectral structure might also be pro- 
duced by dynamical coupling between hydrogen at ad- 
jacent sites. A para-para sensitive tool such as Raman 
or IR may shed further light onto the origin of various 
structures. 




Energy [meV] 

FIG. 3: Our assignments for the observed energy loss peaks. 
The circles and diamonds are the 4 K inelastic neutron scat- 
tering data points of Ref. with the ordinate origins shifted 
arbitrarily. We put cubic splines through the data points as 
a guide to the eye. 
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TABLE II: Calculated translational frequencies (meV/h) of 
adsorbed H2. 





# Ha/cell 


nl 


n2 


n3 


Zn site 


1 


7.1 


15.2 


17.5 


12 


9.2 


12.6 


22.9 


site 


12 


9.5 


14.5 


21.7 




FIG. 4: The adsorption sites. Left panel: our calculational 
results for H2 or D2 without consideration of zero point mo- 
tions. Right panel: the neutron diffraction results [3] for D2 
at 4 K. The dihydrogens shown are located at eight different 
distances perpendicular to the plane of the figure within each 
primitive cell. 



With the rotational and translational energy levels cal- 
culated, we are able to obtain the zero point energy cor- 
rections to the total binding strength, which are about 10 
and 20 meV respectively. Therefore, the effective binding 
strength at low loading (Zn site) is ^100 meV. It is some- 
what larger than the experimental value of 91 meV [2| 
(88 meV in Ref. [H). The slight overestimation of the H2 
binding energy by vdW-DF is also observed for another 
MOF 181. For the O site, our calculated binding includ- 
ing zero point correction is ^90 meV. This latter value 
would probably be smaller relative to the prediction for 
the Zn site, if it were not for a known inaccuracy of an 



extremal limit in the exchange part of the functional [16| . 
The experimental value [2] is ~50 meV. 

Finally, we discuss the hydrogen uptake in this com- 
pound. As shown in Ref. [H, a total of 24 D2/cell can be 
loaded into the system to form a structure as shown in 
the right panel of Fig. [4] where each hydrogen molecule 
is represented by a single point. To make a comparison, 
we put 24 H2 per cell into the MOF and let them re- 
lax. To speed up the relaxation, the dihydrogens were 
initially positioned to satisfy the 3-fold rotation plus in- 
version symmetry, with 6 near each Zn site, 6 near each 
O site, 6 around 3 A from each aromatic ring, and 6 at 
random in the central region. The final result is shown 
in the left panel of Fig. [4] where both hydrogen atoms in 
H2 are displayed. Overall agreement is obtained between 
theory and experiment with the the sorbate-sorbent dis- 
tances in our calculations deviating from experiment by 
about 15%. 

In summary, we studied H2 adsorbed within the 
nanopores of a complex metal organic framework ma- 
terial (MOF- 74) [3] with vdW-DF calculations. We con- 
firmed the primary and secondary binding sites associ- 
ated with unsaturated metal atoms and with triangles of 
O atoms, respectively. The predicted (zero point cor- 
rected) binding energy of the primary site is in good 
agreement with heat of adsorption measurements @,[H at 
low H2 loading. We also identify a third and fourth type 
of site, which are fully occupied at an H2 concentration 
of 24 H2 per primitive cell, in agreement with experi- 
ment [2]. Rotational motion of the adsorbed H2 at the 
primary and secondary sites is analyzed. The dihydro- 
gen behaves as three-dimensional quantum rotor at both 
sites and the hindered rotation splits the triplet state 
into three nondegenerate levels. Based on the calculated 
para-ortho transitions from singlet to triplet, we are able 
to identify the main features in the INS spectra Q for 
different H2 concentrations, with agreement on the 20% 
(5%) level for the primary (secondary) sites. 
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